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We investigate tlie Kondo Lattice Model on 2D clusters using the Finite Temperature Lanczos 
Method. The temperature dependence of thermodynamic and correlations functions are systemati- 
cally studied for various Kondo couplings Jk- The ground state value of the total local moment is 
presented as well. Finally, the phase diagrams of the finite clusters are constructed for periodic and 
open boundary conditions. For the two boundary conditions, two different regimes are found for 
small Jx/t, depending on the distribution of non-interacting conduction electron states. If there are 
states within Jk around the Fermi level, two energy scales, linear and quadratic in Jk, exist. The 
former is associated with the onsite screening and the latter with the RKKY interaction. If there are 
no states within Jk around the Fermi level, the only energy scale is that of the RKKY interaction. 
Our results imply that the form of the electron density of states (DOS) plays an important role in 
the competition between the Kondo screening and the RKKY interaction. The former is stronger if 
the DOS is larger around the Fermi level, while the latter is less sensitive to the form of the DOS. 

PACS numbers: 62.30.-|-d,65.40.-b,66.35.-f a 



I. INTRODUCTION 

The Kondo lattice model (KLM) is used to describe 
compounds containing localized magnetic moments, such 
as heavy fermion systems and Kondo insulatorsi The 
KLM Hamiltonian is given by: 

HkL = 4<TCi,a + Jk'^Ti- Si. (1) 

Here, and are itinerant and local (/-) spins on site 
J, respectively. The KLM takes into account hopping of 
conduction electrons, t, and their Kondo interaction, Jk, 
vi^ith local /-spins. The Kondo term causes a screening 
of local spins, but also induces indirect RKKY interac- 
tions between the local spins on the lattice. These two 
interactions compete, leading to either magnetically or- 
dered or non- magnetic ground states, separated by the 
quantum critical point (QCP).^- The nature of the QCP 
in heavy fermion compounds is one of the central topics 
of condensed matter physics todayi^i^i^iSiL2i2ii£ 

Various analytical and numerical methods have been 
used to study the KLM, like the mean-field approach^^ 
and Quantum Monte Carlo (QMC) method>i2i The finite 
temperature Lanczos method (FTLM)i^ has been used 
to study thermodynamic functions Also the periodic 
Anderson model on finite size clusters has been studied 
with similar methods recently^ 

Here, we use the FTLM and mainly focus the analysis 
on the temperature dependence and the ground state val- 
ues of static correlation functions, of central importance 
for the QCP. We compare onsite vs. intersite correla- 
tions on the lattice and also compare the former with 
the impurity case. The latter is associated with RKKY 
interactions. We also investigate the evolution of the 
total onsite moment with Kondo coupling strength Jk. 
The main advantage of the FTLM is that it involves no 
uncontrolled approximations and is exact for finite size 




FIG. 1: DOS for the tight-binding model in 2D is shown 
shaded. Unshaded bars represent discrete levels of the 2D 8- 
site cluster (for Jk = 0) with PBC. The black bars represent 
levels with OBC. The bar heights correspond to the level de- 
generacy (divided by 8). Left Inset: The 8-site cluster with 
site indices. Right Inset: The first Brillouin zone with the 
k-points corresponding to Bloch states of 8-site cluster. The 
shaded part represents the Fermi volume for half filled case. 
There are 6 degenerate states on the "Fermi surface" when 
the multiplicities are taken into account correctly. 



clusters. It treats both, Kondo screening and RKKY in- 
teractions, correctly. 

We consider the 8-site cluster of a square lattice, as 
shown in the inset of Fig. ^ with periodic (PBC) and 
open (OBC) boundary conditions at half filling. 

We will show that for this cluster the two different 
boundary conditions lead to two interesting physical sit- 
uations, which may have relevance to the two types of 
QCP's observed (cL Ref.Qand references cited therein). 
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There is of course only one solution of the KLM in the 
thermodynamic limit. However, the real materials have 
more complex band structures than described by the n. n. 
hopping in (^). By varying boundary conditions for the 
KLM on the small cluster one may mimic band structures 
with different characteristic behaviour of the density of 
states (DOS) around the Fermi level. For this reason we 
present in detail the analysis of the KLM on the finite 
cluster for the PBC and the OBC and compare the two 
cases. In Sec. II we discuss the specific heat anomalies 
obtained in the two cases. Sec. Ill discusses the central 
problem of Kondo-lattice physics: The competition and 
crossover behaviour of onsite and intersite correlations 
as function of temperature and control parameter Jk/^. 
The associated Kondo screening behaviour evident from 
the susceptibility and total local moment is presented in 
Sec. IV. Finally Sec. V summarizes the results in a phase 
diagram and gives the conclusions. 



II. SPECIFIC HEAT 

The specific heat for the PBC and the OBC, for various 
coupling constants Jk, is shown in Fig. [21 

The low-temperature specific heat for non-interacting 
case, Jk = 0, with PBC deviates significantly from the 
thermodynamic limit. The reason is that 6 out of 8 elec- 
trons are placed on the " Fermi surface" , shown in the 
inset of Fig. ^ thus the ground state is highly degen- 
erate. The specific heat for OBC is closer to the ther- 
modynamic limit. The peak from the hopping term for 
Jk = is much more pronounced and at lower tempera- 
ture, rO^^^o, for OBC than for PBC. The electronic states 
for OBC are no longer Bloch states, but the important 
point is that the electron energies are shifted from, and 
distributed below, the " Fermi level" . The distribution of 
these discrete levels is such that it relatively well mim- 
ics the continuum DOS in the thermodynamic limit, as 
shown in Fig. ^ 

The specific heat for PBC and OBC, when Kondo in- 
teraction is turned on, shows qualitatively different be- 
havior in the weak coupling regime, Jk/^ ^ 2, as shown 
in Fig. 121 For PBC, a double peak structure is observed 
at low temperatures for 0.5 < Jk/^ ^ 2. We may denote 
the positions of the upper and lower peak by T^aiFi ^^'^ 
^ma3?2; respectively. For OBC only one low temperature 
peak appears at Tma^S.- The origin of these peaks are 
spin degrees of freedom, whose degeneracy is lifted by 
Kondo interaction. In the strong coupling limit there is 
only one peak in the specific heat, and it does not depend 
on the boundary conditions. 

The specific heat for OBC agrees with the results ob- 
tained with the same method for the 10-site cluster with 
PBC inRef.Q as well as with the results obtained with 
the QMC in Ref.[li 
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FIG. 2: (Color online) The specific heat calculated for various 
values of Jx/i, as indicated in the legend. The legend is 
valid for all the figures in this work, if not explicitly specified. 
Above: for PBC. The two low temperature peaks can be seen 
for 0.5 <JK/t< 2. In the inset the lower peaks are magnified. 
Below: for OBC. The JK/t~ curve agrees well with the 
thermodynamic limit of the Hamiltonian Q. For JK/t> 
there is only one low temperature peak. 



III. SPIN CORRELATION FUNCTIONS: 
ONSITE SINGLET FORMATION VS. RKKY 
INTERACTIONS 

The nature of the specific heat peaks may be revealed 
by considering appropriate static correlation functions. 
We focus the analysis on the following two: (ti ■ Si) is 
the onsite correlation function between local and itiner- 
ant spins, while (Si • S5) is the intersite correlation func- 
tion between local spins on neighboring sites for PBC. 
For OBC these are {t^ ■ S3) and (S3 • Sg), i.e., the in- 
ner sites 3 and 6 are considered, because they are more 
representative than sites on the edges of the cluster (cf. 

Fig. mi. 

A. Periodic Boundary Conditions 

The onsite correlations, (ri • Si), for PBC as function 
of temperature are shown in Fig. (21 They are also com- 
pared with the onsite correlations for the impurity case. 
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FIG. 3; (Color online) The onsite correlations in the lat- 
tice (full lines) as function of temperature and the compar- 
ison with the impurity case (broken lines), for jR/i values 
as indicated in the legend. The correlation strength increases 
monotonically with increasing Jx/i- In the impurity case, the 
impurity local spin is placed on site 1 for the PBC, and site 
3 for the OBC. Above; for PBC. Below: for OBC. Due to 
the small system size for the impurity case, there are some 
unphysical deviations at low temperatures and weak coupling 
JkA for the OBC. 
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FIG. 4: (Color online) Correlation functions for PBC. Above: 
Normahzed (ti - Si) as function of temperature divided by 
Jk. Below: Normalized (Si • S5) as function of temperature 
divided by the square of Jk. The normalization is achieved 
by dividing the curves with the ground state values. In the 
inset the ground state (T = K) correlations as function of 
Jx/t are shown. 



where there is only one local spin, at site 1, on the 8-site 
cluster. In the weak coupling regime the onsite screening 
is decreased in the lattice, because there are less itinerant 
spins per local spin available for the screening. 

The inflection points of the correlation functions, as 
function of temperature, indicate that the upper peak in 
the specific heat, at T^f^^i, corresponds to the forma- 
tion of the onsite correlations, whereas the lower one, at 
^ma3?2i corresponds to the formation of the intersite cor- 
relations. The ground state {T — K) correlations for 
PBC as function of coupling strength Jx/i are shown in 
the inset of Fig. ^ It is seen that the onsite correlations 
are stronger than the intersite correlations for all Jk/^- 
Another characteristic is that (xi • Si) jumps to the large 
absolute value as soon as Jk is turned on. This is a con- 
sequence of the degenerate itinerant degrees of freedom 
at Fermi level for Jk = 0. The free itinerant spins on 
the " Fermi level" screen the local spins onsite as soon as 
Jk > 0; the ground state splits off and is characterized 
by finite (ti - Si) correlations. The unscreened parts of 



the local spins form intersite correlations, (Si • S5), below 
^ma3?2; ^uc to thc RKKY interaction. The magnitude of 
(Si • S5) increases rapidly as Jk is turned on. Because 
of the very small energy scale (quadratic in Jk), conver- 
gence problems appear in the calculation of ground state 
(Si • S5) correlations for Jk/^ < 0.5 and we were not 
able to determine whether they also jump for finite Jk- 
In the strong coupling regime the intersite correlations 
vanish, while thc onsite correlations saturate at —3/4, as 
the onsite singlets arc formed. 

The correlation functions, normalized to the ground 
state values, versus scaled temperature, as done in 
Ref. fl^ are shown in Fig. 0] The onsite correlations 
are characterized by a single temperature scale, perfectly 
linear in Jk/^. The intersite correlations are also char- 
acterized by a single temperature scale proportional to 
Jk^ /t. This shows that for PBC the onsite screening 
is independent of the intersite correlations, which may 
be associated with RKKY interactions. This is a conse- 
quence of the special symmetry of the 8-site cluster which 
implies a large fraction of conduction electron sitting on 
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FIG. 5: (Color online) Correlation functions for OBC. Above: 
Normalized (ts -83) as function of temperature divided by 
Jk- Below: Normalized (S3 ■ Se) as function of temperature 
divided by Jy^ jt. The normalization is achieved by dividing 
the curves with the ground state values. In the inset the 
ground state (T = K) correlations as function of Jx/t are 
shown. 



the " Fermi surface" . Nevertheless it is an interesting re- 
sult which may have relevance for a case where the con- 
duction electron DOS is strongly peaked at the Fermi 
level. 



B. Open Boundary Conditions 

The onsite correlations, (X3 • S3), for OBC as function 
of temperature and the comparison to the impurity case 
are also shown in Fig. |3| In contrast to the PBC, for 
OBC the onsite correlations are enhanced in the lattice. 
In the case of OBC there are no "free" itinerant spins, i.e., 
conduction electrons on the " Fermi level" . In the weak 
coupling limit the onsite screening is not a direct process 
and is weak. The interesting point is the appearance of a 
new low energy scale in the lattice, which may be noted 
as a low-temperature downturn in the onsite correlations 
in Fig.El 

The ground state correlations are shown in the in- 
set of Fig. O The onsite correlations, (T3 ■ S3), do not 
jump as Jk is turned on, but rather increase continuously 



with Jk. The intersite correlations, (S3 • Sg), are much 
stronger than (Si ■ S5) for PBC and in the weak couphng 
limit they are even larger in absolute values than the on- 
site (t3 • S3) correlations. It is understandable, because 
the local spins are only weakly screened, due to the lack 
of "free" itinerant spins. Therefore, the unscreened local 
spins can form strong intersite correlations. 

The normalized correlations as function of the scaled 
temperature for OBC are shown in Fig. |S1 The (r3 • S3) 
scale with Jk only in the strong coupling limit. In the 
weak coupling limit two characteristic temperatures or 
energy scales can be identified, as already mentioned 
above. The upper one is associated with the energy split- 
ting of the discrete energy levels (cf. Fig.^, and is clearly 
a finite size effect. The lower energy scale is associated 
with the additional onsite screening, due to the presence 
of other local spins in the lattice. The intersite (S3 • Sg) 
correlations are still approximately governed by a single 
energy scale, Jk^/^, just as for PBC. Indeed the behavior 
of the normalized intersite correlation functions for OBC 
and PBC are quite similar, despite of the fact that the 
magnitude differs by a factor of two. 

IV. MAGNETIC SUSCEPTIBILITY AND 
TOTAL LOCAL MAGNETIC MOMENT 

The spin correlations are also directly reflected in mag- 
netic susceptibility, shown in Fig.|Bl for PBC and OBC. It 
is multiplied by Jk and plotted versus temperature scaled 
with Jk, as was done in Ref.fl^ The two peaks for PBC 
are clearly seen for small Jk/^, as in the specific heat. 
However, in the magnetic susceptibility even the lower 
peak for JK/i= 0.5 is clearly seen. This is because the 
lower peaks in the magnetic susceptibility are much more 
pronounced and at slightly higher temperatures than the 
corresponding peaks in the specific heat. The position 
of the lower peak scales well with J^ (not shown in the 
figure). The lower peak is not distinguishable any more 
for Jk/* ^ 2. The position of the higher peak scales with 
Jk. For OBC only one large peak is seen. For Jk/^ ^ 2 
it scales approximately with J"^. In the strong coupling 
limit, Jk/^ ^ 2, the peak scales linearly with Jk, just 
as for PBC. The results for OBC are similar to the re- 
sults for the 10-site cluster with PBC in Ref. Q These 
results are also similar to the results of the QMC calcu- 
lations for larger clusters in Ref. However, in Ref.H^ 
the magnetic susceptibility peak scales with Jk^, while 
in the onsite correlations the characteristic temperature 
is proportional to Jk, suggesting the existence of two low 
energy scales. 

In the insets of Fig. the ground state values of the 
total local moments are shown, defined as: 

(mL) = ((r, + S,)^) = (rf) + (S?) + 2(r, ■ S,). (2) 

The interesting quantity here is (^t^) which is a measure 
of itinerant spin localization. In the strong coupling limit 
(rf) — 3/4, i.e., the itinerant spins become localized as 
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FIG. 6: (Color online) The scaled magnetic susceptibility, 
as indicated on the axes. Above: for PBC. Below: for OBC. 
In the insets the total local moments in the ground state (at 
T = K) as function of Jx/i are shown, (/xfoc) ~ {(''"» + Si)^), 
(i = 1,3 for the PBC and OBC, respectively). 



the onsite singlets are formed and the hopping of conduc- 
tion electrons is frozen. When the conduction electrons 
are completely delocalized, i.e., for Jk= 0, {t^) = 3/8, 
because the probability of finding one electron with spin 
up or down on one particular site is 1/2. Accordingly, 
(m?oc) = (m?oc)o = 9/8 for Jk = while for Jk oo 
it vanishes, as the onsite singlets are formed and the lo- 
cal moments completely screened. There is a jump for 
JK/t> for PBC. It has the same origin as the jump for 
the onsite correlations (ri-Si), i.e., the highly degenerate 
unperturbed conduction electron ground state. A jump 
is also present in (t^). For OBC we see a continuous 
decrease of {Hi^c)- 

The (/^i^oc) ^ good quantity to define a critical value 
of Jk/^ at which the crossover from predominantly mag- 
netic to predominantly nonmagnetic ground state occurs. 
The other possibility would be to compare directly the 
intersite and onsite correlations, but the intersite cor- 
relations depend strongly on the cluster size. There- 
fore we define a critical (Jk/^)c as the value for which 
{<c) - (a*1oc)o/2. Then for OBC, (JK/t)c « 1-4 (cf. 
Fig. , in good agreement with the values obtained from 
the Quantum Monte Carlc^-. and the mean-field analy- 
sisM. For PBC it would be 0, but for this small cluster 



FIG. 7: (Color online) The phase diagrams for the 8-site 2D 
clusters. In the insets the weak coupling region is magnified. 
Above: for PBC. The inflection points are approximately de- 
termined as the temperatures where the normalized (Si • S5) 
has the value of -0.75. The dashed line is the fitting curve 
assuming quadratic dependence on Jk. The inflection points 
of the onsite correlations agree with the positions of the upper 
peak in the specific heat and are not shown in the Figure. Be- 
low: for OBC. The inflection points are approximately deter- 
mined as the temperatures where the normalized correlation 
functions achieve the values of -0.8 and -0.75 for (ts • S3) and 
(S3 ■ Ss), respectively. 



size, OBC are much more representative for the thermo- 
dynamic limit of the KLM. 



V. PHASE DIAGRAM AND CONCLUSIONS 

We may summarize the results in form of phase dia- 
grams for PBC and OBC, shown in Fig. [7| In the phase 
diagrams the positions of specific heat maxima are plot- 
ted as function of coupling strength, Jk/^. Also the in- 
flection points of the correlation functions are plotted, 
and they mostly follow the positions of the correspond- 
ing specific heat maxima. 

In the weak coupling regime, for PBC, two character- 
istic temperatures, T^f^i and are found. The 
^ma3?2 cannot be determined for Jk/*^ 2, but we make 
the continuation with the (Si • S5) infiection point tem- 
perature, T^Si-Ss); which can be defined for all coupling 
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strengths. Therefore we conclude that the correspond- 
ing anomaly in the specific heat is still there, but is too 
small to be identified. The inflection points of the cor- 
relation functions reveal that the characteristic temper- 
atures T^^x,i s-iid ^max^2 corrcspond to the energy scales 
of onsite correlations and intersite correlations, respec- 
tively. The former is linear in Jk for all J^/t. The latter 
is proportional to Jy^ jt^ and may be associated with the 
RKKY interaction. 

For OBC, there is only one characteristic temperature, 
^max"i: ^ Jk/^- In the weak couphng regime, for 
JY,/t< 2, below T^^^i predominantly intersite correla- 
tions are formed. However, T^^-^i is also the charac- 
teristic temperature of onsite correlations. This is seen 
from the positions of inflection points for the correspond- 
ing correlations, also shown in the lower part of Fig. [7| 
Because the intersite correlations dominate and the char- 
acteristic temperature is approximately quadratic in Jk, 
the corresponding energy scale, in the weak coupling 
regime, may be associated with the RKKY interaction. 

It is important to note that T'^^i m T^^^2j the 
RKKY interaction is not very sensitive to the distribu- 
tion of non-interacting conduction electron states. This is 
plausible with the known expression for the RKKY char- 
acteristic temperature, Trkky oc J^/W (e.g.^), where W 
is the band width. It is also important that the RKKY 
energy scale is the energy scale of the onsite correlations 
in the weak coupling regime for the OBC. We may gen- 
erally conclude that the RKKY energy scale is the low- 
est energy scale of the KLM; if the local spins are not 
screened above Trkky7 then it is the only energy scale 
of the system. 

In contrast, the onsite correlations, i.e., the onsite 
screening of local spins, is very sensitive to the number of 
states within Jk around the " Fermi level" , as seen from 
the difference between T^^^^^ and T^^^^^ for PBC and 
OBC, respectively. This is plausible with the expres- 
sion for the single impurity Kondo temperature, which 
depends exponentially on the DOS at the Fermi level. 



p{<^f); Tk = epexp [-l/(Jp(eF))]. 

In the strong coupling regime, only one characteristic 
temperature is found for both PBC and OBC. It cor- 
responds to the energy of the onsite singlet formation, 
which is linear in Jk . Boundary conditions are unimpor- 
tant, because the local physics determines the properties 
of the system. 

We may argue on the relevance of these results for the 
system in the thermodynamic limit. Recently a QCP was 
found5*S*& whose properties deviate from the expected be- 
havior (e.gAiS). The properties of the system at the 
QCP are determined by the competition between the on- 
site screening and the RKKY interaction. Here we have 
shown that the competition between the onsite screening 
and the RKKY interaction strongly depends on the par- 
ticular form of the distribution of non-interecting con- 
duction electron levels. We may anticipate that in the 
thermodynamic limit this means a strong dependence on 
the particular form of non-interacting conduction elec- 
tron DOS. A strongly peaked DOS near the Fermi level 
favors onsite screening, while the RKKY interaction is 
not very sensitive to the particular shape of the DOS. 
These results indicate that the shape of the DOS is very 
important for the properties of the heavy fermion systems 
in the vicinity of the QCP. 

Finally, the qualitative behavior of the correlation 
functions for the KLM with OBC is similar to the KLM 
without charge degrees of freedom,^'' which is a general- 
ized Kondo necklace model in 2D from Doniach's original 
work.- This shows that this simplified model describes 
well the spin degrees of freedom in the KLM at half- 
filling. 
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